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THE REST-FRAME FUV MORPHOLOGIES OF STAR-FORMING GALAXIES AT Z - 1.5 AND Z - 4 
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ABSTRACT 

We apply a new approach to quantifying galaxy morphology and identifying galaxy mergers to 
the rest-frame far-ultraviolet images of 82 z 4 Lyman break galaxies (LBGs) and 55 1.2 < z < 1.8 
emission-line galaxies in the GOODS and Ultra Deep Fields. We compare the distributions of the Gini 
coefhcient (G), second-order moment of the brightest 20% of galaxy light (M20), and concentration (C) 
for high-redshift and low-redshift galaxies with average signal to noise per pixel > 2.5 and Petrosian 
radii >0.3". Ten of the 82 LBGs have M2Q > —1.1 and possess bright double or multiple nuclei, 
implying a major-merger fraction of star-forming galaxies 10-25% at Mpuv < ~20, depending on 
our incompleteness corrections. Galaxies with bulge-like morphologies (G > 0.55, M20 < —1.6) make 
up ~ 30% of the z ~ 4 LBG sample, while the remaining 50% have G and M20 values higher than 
expected for smooth bulges and disks and may be star-forming disks, minor-mergers or post-mergers. 
The star- forming z ~ 1.5 galaxy sample has a morphological distribution which is similar to the UDF 
z ^ 4 LBGs, with an identical fraction of major-merger candidates but fewer spheroids. The observed 
morphological distributions are roughly consistent with current hierarchical model predictions for the 
major-merger rates and minor-merger induced starbursts at z ~ 1.5 and ~ 4. We also examine 
the rest-frame FUV — NUV and FUV — B colors as a function of morphology and find no strong 
correlations at either epoch. 

Subject headings: galaxies:evolution - galaxies: high-redshift - galaxies:interacting ~ galaxies:structure 
- ultraviolet: galaxies 



1. INTRODUCTION 

Galaxy interactions and mergers are expected to play 
a dominant role in the formation of stars and evolu- 
tion of galaxies, particularly in the early universe. The 
present hierarchical models of galaxy assembly and star 
formation fail to predict enough red spheroidal galax- 
ies - the assumed end products of major galaxy mergers 
(e.g. Toomre 1977; Mihos & Hernquist 1996) - at both 
early and late times (e.g. Cimatti et al. 2004). Observa- 
tions of on-going mergers at high redshift and the earliest 
spheroidal systems are crucial to constraining the role of 
mergers in the evolution of galaxies and the build-up of 
stellar mass. One of the most accessible ways of identi- 
fying galaxy mergers and their end-products is by their 
morphologies. 

Morphological studies of very distant objects have be- 
come possible in the last decade with the high-resolution 
capabilities of the Hubble Space Telescope. Galaxies 
at redshifts greater than one appear to be systemati- 
cally more disturbed than local galaxies (Abraham et 
al. 1996), implying strong evolution in the number den- 
sity of merging and interacting systems (Abraham & 
van den Bergh 2002). The first HST images of galax- 
ies above z ~ 3 identified via the Lyman break tech- 
nique (Steidel et al. 1995; Madau et al. 1996) showed 
a ranges of morphologies, from compact proto-bulges to 
disturbed major-merger candidates (Giavalisco, Steidel, 
& Macchetto 1996; Lowenthal et al. 1997; Dickinson 
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1999; Giavalisco 2002). In recent years, high redshift ob- 
jects have been identified by other techniques and subse- 
quently studied with HST. These include the extremely 
red objects (EROs) and sub-mm sources at z ~ 1 — 3. 
EROs possess a range of morphologies from spheroidal 
to irregular (Moustakas et al. 2004), with an increasing 
fraction with multiple nuclei at z > 2 (Daddi et al. 2003). 
Optical HST images of a small sample of sub-mm sources 
show that 70% possess multiple components (Chapman 
et al. 2003). 

The connection between the morphologies of very early 
objects and those observed at later times is not straight- 
forward. In the current picture of galaxy evolution, 
a galaxy's morphology is not constant through cosmic 
time, nor does it evolve along a single linear path. Star- 
forming disks are generally grown through the accretion 
of cold intergalactic gas, while bulges and spheroids are 
typically produced in mergers and interactions. For ex- 
ample, in one simulation of the evolution of a present-day 
elliptical, its progenitors start as disks grown at redshifts 
> 4 (Steinmetz & Navarro 2002). These primordial disk 
galaxies merge to form a bulge at z ^ 3. The merger 
remnant then dynamically relaxes, accretes more gas and 
reforms a disk between l<z<2. Atz<l, the galaxy 
undergos yet more interactions and minor-mergers until 
its gas supply is exhausted and it resembles a red ellip- 
tical at z ~ 0. 

Despite the complex evolutionary path a single object 
may take, the evolution of the morphological distribu- 
tion of galaxies and the number of merging systems at 
z > 1 is an important tracer of the formation of the 
Hubble sequence. While disks are quite fragile and can 
be destroyed and then regrown, numerical simulations 
suggest that, once created, bulges are not transient and 
only grow in mass as a galaxy evolves (Mihos & Hern- 
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quist 1996, Stcinmctz & Navarro 2002). At z > 3, the 
majority of galaxies are predicted to be disk-dominated 
(Somerville, Primack, & Faber 2001, hereafter SPFOl). 
However, many of these proto-disks may not resemble 
local spirals because the disks might not be stable (Mao, 
Mo, & White 1998) and because a substantial fraction 
are expected to be undergoing minor mergers (SPFOl). 
Also, by z 3, approximately 20% of objects may 
have already transformed into bulge-dominated systems 
(SPFOl). The fraction of dark matter halos undergoing a 
major merger is expected to increase as a power-law pro- 
portional to (1 + z)'' out to at least z ^ 2. Above z ^ 2, 
theoretical predictions imply that the merger rate flat- 
tens such that ^ 20—30% of massive dark matter lialos at 
z> 2 should merge within a Gyr (Gottlober, Klypin, & 
Kravtsov 2001). High resolution hydrodynamical simula- 
tions indicate that the associated gas-rich galaxy mergers 
create bulge-dominated galaxies (e.g. Mihos & Hernquist 
1996). 

Some of these mergers shoifld evolve into the z < 1 red 
spheroidal population, particularly if they live in dense 
environments. The spatial clustering of high redshift ob- 
jects has been used to argue that bright Lyman break 
galaxies are the likely progenitors of highly-clustered 
present-day ellipticals (Giavalisco et al. 1994; Giavalisco 
2002), although fainter LBGs are much less clustered 
(Giavalisco & Dickinson 2001) and alternative models 
are possible (e.g. Wechsler et al. 2001). Other z > 2 
galaxy populations are also highly clustered (e.g. Daddi 
et al. 2001). Yet it remains ambiguous to what extent 
the observed high redshift galaxies track the progenitors 
of local spheroids, since their stellar masses, mctallici- 
ties, and star-formation rates must evolve substantially 
at later times (Dickinson et al. 2003; Ferguson et al. 
2002; Lowenthal et al. 1997). 

Despite the great advances in identifying high-redshift 
galaxies and obtaining high-resolution images, the dis- 
tribution of intrinsic galaxy morphologies at early times 
is still uncertain. Complicating the issue is the wave- 
length dependence on observed morphology. The high- 
est spatial resolution HST images are those taken in op- 
tical wavelengths with the Advanced Camera for Sur- 
veys (ACS). At z > 1.2, ACS can only sample rest-frame 
wavelengths less than the 4000A break, where dust and 
young star-forming regions dominate the morphologies 
of local galaxies (Giavalisco et al. 1996; Hibbard and 
Vacca 1997; Goldader et al. 2002). Therefore the mor- 
phological distributions of high-redshift objects classified 
by their rest-frame ultraviolet morphologies should not 
be compared to the rest-frame optical morphological dis- 
tributions of nearby objects. 

A second problem is the intrinsic difficulty of mea- 
suring the morphologies of faint and compact objects 
like Lyman break galaxies. Lyman break galaxies at 
z 4 have typical rest-frame far-ultraviolet (FUV) half- 
light radii between 0.1 and 1.0 " (Ferguson et al. 2004; 
Bouwens et al. 2004), therefore many LBGs with sizes 
below at few tenths of an arcsecond are not spatially 
resolved enough to measure their morphologies even in 
ACS images. The compact sizes may bias standard 
morphological analyses that are commonly applied to 
lower redshift objects, such as concentration measures 
and bulge-disk decompositions. Also, images of faint 



high-redshift objects are likely to miss light in the outer, 
lower-surface brightness regions of the objects. This 
missed light will bias the observed sizes low, and also af- 
fect the measured concentrations and surface-bright pro- 
files. Finally, morphology measures like rotational asym- 
metry and dumpiness which involve subtracting images 
and/or a correction for the background noise signal are 
less robust at low signal-to-noise levels (Lotz, Primack, 
& Madau 2004). 

A third possible problem is the assumption that dis- 
turbed morphologies at high redshift directly correlate 
with mergers and interactions. Because of the roughly 
500 million year timescales associated with merger ac- 
tivity, we do not truly know if a disturbed object has 
undergone a recent merger from a snapshot of its mor- 
phology alone. The observed morphology depends on 
the observational effects (rest-frame wavelength, view- 
ing angle, dust extinction) and the physical conditions of 
the merger/interaction (mass ratio, orbital parameters, 
gas fractions, bulge fractions, internal kinematics). Spec- 
troscopy can sometimes confirm the presence of two or- 
biting galaxies (e.g. the Superantenna, Mihos & Bothim 
1998) or tidal in-fall (e.g. IRAS 15204-3342, Arriaba & 
Colina 2002). Recent observations of morphological ir- 
regular, elongated systems at z 2 have failed to show 
the significant rotation that is expected at early merger 
stages (Erb et al. 2004), and recent studies of kinematic 
close pairs at z ^ 1 indicate a lower merger rate than 
that derived using morphological asymmetry (Lin et al. 
2004). However, the detection of kinematic merger signa- 
tures will depend on the merger stage as well as viewing 
angle and the slit-position. ULlRGs in the final merger 
stage can exhibit complex velocity fields without obvi- 
ous rotational centers (eg. IRAS 23128-5919, Mihos & 
Bothun 1998). While this paper assumes that morpho- 
logically identified double- nucleated galaxies are merging 
systems, theoretical and empirical calibrations of quanti- 
tative morphologies and kinematic signatures as merger 
indicators are important tasks for future work. 

In this paper, we apply a new approach to quantifying 
galaxy morphology and identifying mergers at high red- 
shift. First presented in Lotz, Primack and Madau (2004; 
hereafter LPM04), this technique measures the relative 
distribution of galaxy pixel flux values or Gini coefficient 
(Abraham et al. 2003) and the second-order moment of 
the brightest 20% of the galaxy's pixels (M20). These 
two quantities are extremely robust to decreasing signal- 
to-noise and change less than 10% for average signal- 
to-noise per pixel > 2 at spatial resolutions better than 
500 pc per resolution element (LPM04). At rest-frame 
optical wavelengths, we found that local galaxy mergers 
separate cleanly from the local galaxy Hubble types in 
Gini coeflttcicnt (G) and M20. 

Here we present the non-parametric rest-frame FUV 
morphologies for 82 z ~ 4 Lyman break galaxies ob- 
served in the GOODS and UDF survey and compare the 
morphology distributions to a spectroscopically selected 
sample of 55 z ~ 1.5 galaxies in the GOODS fields. In 
§2, we discuss the selection of both galaxy samples and 
note that our analysis is limited to FUV bright star- 
forming galaxies at both epochs. In §3 we describe our 
non-parametric measures, and how these are applied to 
the GOODS and UDF images. The size and surface- 
brightness selection effects are examined in §4 with both 
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a comparison of UDF and GOODS-South morphologies 
and a scries of simulations. The observed morphological 
distributions of objects with sufficient sizes and signal- 
to-noise are presented in §5. We also compare the galaxy 
morphologies to their colors for a sub-set of objects in the 
GOODS-South and UDF. Finally, in §6 we discuss the 
morphological distributions in the context of hierarchical 
model predictions for star- forming galaxies at z > 1.5. 

2. GALAXY SAMPLES 

Images with high spatial resolutions (better than 1 kpc 
per resolution element) are needed to reliably measure 
galaxy morphologies. Assuming the concordance cosmo- 
logical model {Ho = 70 km s~-^ Mpc^^, = 0.3, = 
0.7), galaxies at rcdshifts greater than z ^ 1 have angu- 
lar scales ^ 5 kpc per arcsec. The point spread hmction 
(PSF) of the Advanced Camera for Surveys (ACS) on the 
Hubble Space Telescope {HST) has a FWHM = 0.12 ", 
or > 800 pc at 2; > 1. Two recent ACS programs 
have provided an archive of deep, high-resolution, multi- 
wavelength images that are invaluable for the study of 
very distant galaxies: Great Observatories Origins Deep 
Survey (GOODS; Giavalisco et al. 2004) and the Ul- 
tra Deep Field (PI: S. Beckwith). The GOODS ACS 
images cover two 160 □' fields (the Hubble Deep Field 
North and the Chandra Deep Field South) in four band- 
passes (F435W, F606W, F775W, and F850W), with a 10 
(7 detection hmit of ~ 26.7 in F850W. The UDF covers 
a much smaller area (~ 11.3 □') in the Chandra Deep 
Field South, but is observed in the same ACS band- 
passes and has a deeper limiting magnitude 28.4 in 
F850W). An additional advantage of these two surveys 
is that the overlap between the extremely deep UDF and 
somewhat shallower GOODS fields allows us to empiri- 
cally determine the limiting signal-to-noise for our mor- 
phology measures. 

We use these two sets of ACS images to measure the 
morphologies of Lyman break galaxies (LBGs) a,t z ^ A 
and a lower-redshift comparison sample at 2; ~ 1.5. At 
z ^ 4, even the longest wavelength z-band (with a cen- 
tral wavelength Ac = 8500A) ACS images probe light 
emitted at rest-frame far-ultraviolet (FUV) wavelengths. 
Longer wavelength NICMOS images do exist for the UDF 
and the HDFN, but these have significantly worse reso- 
lution (NIC3 PSF FWHM > 0.22 - 0.35 ")• We wiU 
restrict our analysis to the FUV morphologies of the Ly- 
man break galaxies and the comparison sample of lower- 
redshift galaxies at 1.2 < z < 1.8. 

2.1. 2; ~ 4 Lyman-break Galaxies 

The z ^ 4: LBGs were selected from the GOODS 
photometric catalogs using the "dropout" technique (eg. 
Madau et al. 1996, Steidel et al. 1996), via the B435- 
dropout color-criterion of Giavalisco et al. (2004). This 
color-cut is sensitive to star-forming galaxies with low 
to moderate extinctions and rcdshifts 3.4 < z < 4.5 
with an expected median redshift z ^ 3.9. The deepest 
GOODS images are in the V (F606W) and i (F775W) 
bands, which are sufficiently redward of the Lyman break 
at 2 ~ 4 to detect and quantify the LBG morphologies. 
We have chosen to examine the LBG morphologies in 
the summed V and i ACS images, in order to boost the 
signal-to-noise per pixel for each galaxy and obtain the 



largest possible sample of LBGs with measurable mor- 
phologies. The effective wavelength of these summed 
images is 6272 A, or rest-frame 1425 A at 2: = 3.9. We 
have compared the sizes and quantitative morphologies 
of high signal-to-noise LBGs measured in the V + i im- 
ages to those measured in summed i + z images, and 
have found no evidence for patchy Lyman- alpha for- 
est absorption in the V band images which might have 
affected the apparent morphologies. As high spatial 
resolution is not necessary to obtain reliable photome- 
try, we also utilize the European Southern Observatory 
(ESO) VLT/ISAAC ground-based near- infrared photo- 
metric survey of the GOODS-S (CDFS) fields (Vandame 
et al. 2004) to compare the rest-frame FUV morpholo- 
gies to rest-frame FUV — B colors (§5.2). 

Each LBG candidate brighter than 26.6 mag in the 
summed V + i images was visually inspected, and ob- 
jects which were spurious detections (i.e. saturated star 
diffraction spikes, cosmic-rays) were removed from the 
sample. A fraction of the LBG candidates are compact 
and red enough to be late- type M stars. However, as 
discussed below, all objects with < 0.3" were removed 
for our morphological analysis and we expect no stellar 
contamination for our subsample of LBGs. The GOODS 
LBG catalog has ~ 660 B-dropouts detected above 10 a 
in the summed V + i images (^ 26.6 mag). The smaller 
but deeper UDF LBG catalog has 77 B-dropouts de- 
tected above 10 a in the summed V + i image (~ 27.9 
mag). 

2.2. 2 1.5 galaxies 
At 2 ~ 1.2 - 1.8, the bluest GOODS ACS B (F435W) 
images sample the rest-frame FUV at ~ 1980 - 1550 A, 
slightly redder than V + i images for the -B-dropouts. 
Identifying objects in this redshift range is notoriously 
difficult. Spectroscopic rcdshifts are difficult to obtain 
because few emission lines or spectral breaks arc red- 
shifted into the observed optical wavelengths, and opti- 
cal color selection fails because the Lyman break is bluer 
than the optical band-passes. Photometric rcdshifts suf- 
fer from similar problems and we found that the incidence 
of catastrophic failure (|(zspec - Zphot)\/ + Zspec) > 0.2) 
in this redshift range was 25%. Nevertheless, spectro- 
scopic confirmation for some GOODS z ^ 1.5 galaxies 
has been obtained by several surveys. The public Trea- 
sury Keck Redshift Survey with DEIMOS on the Keck II 
telescope has obtained rcdshifts for 1440 objects out of 
a Rab — 24.4 magnitude-limited sample of 2018 galax- 
ies in the GOODS-North field (Wirth et al. 2004). The 
TKRS catalog has "good" quality rcdshifts for 112 galax- 
ies between 1.2 < 2 < 1.8. The VIMOS VLT Deep Sur- 
vey has publicly released rcdshifts for 1599 objects with 
Iab < 24 mag in the GOODS-South (CDFS) field ob- 
tained with VIMOS at the ESO Very Large Telescope 
(Le Fevre et al. 2004). Out of the VIMOS VLT Deep 
Survey catalog objects, 68 galaxies had good quality 
spectra with 1.2 < 2 < 1.8. Finally, the ESO /GOODS 
spectroscopic program with the F0RS2 spectrograph at 
the VLT has obtained spectroscopic rcdshifts for 243 ob- 
jects with 2850 < 24.5 mag in the GOODS-South field 
(Vanzella et al. 2004), and provides an additional 53 ob- 
jects with high quality rcdshifts between 1.2 < 2 < 1.8 to 
our sample. The morphologies for all 408 1.2 < 2 < 1.8 
objects were measured in the GOODS ACS B images. 
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Wc note that all of these spectroscopic surveys are bi- 
ased towards star-forming galaxies with strong emission 
lines. Spectroscopic redshifts are much more difficult to 
obtain for red galaxies at these redshifts, therefore red or 
weak emission-line galaxies are likely to be missing from 
this sample. 

3. MORPHOLOGICAL ANALYSIS 
3.1. Pre-processing 

The GOODS and UDF ACS images are large mosaics 
containing thousands of galaxies. In order to identify the 
galaxies of interest and mask out foreground/background 
objects, wc ran the galaxy detection and photometry 
software SExtractor (Bertin & Arnouts 1996). Segmen- 
tation maps were created for each mosaic image {V + i 
or B) using SExtractor with the detection threshold set 
to 0.6 a and a minimum detection area of 16 pixels. The 
resulting segmentation maps were used as masks for fore- 
ground/background objects. 

Some of the LBGs have very close neighbors. Bright 
star-forming regions of local galaxies may appear as 
several distinct objects when viewed in the PUV (e.g. 
Goldader et al. 2002). To prevent the artificial break-up 
of our LBGs, we compared the SExtractor segmentation 
map for each LBG with the B — V and V — i colormaps. 
If the pixels of a nearby object meeting the LBG color 
criteria overlapped the pixels assigned to the LBG, or if 
another LBG from the GOODS B-dropout catalog was 
within 1.5 ", the segmentation map was edited to include 
that object's pixels with the original LBG. This affects 
4 objects in our final UDF LBG sample and 7 objects in 
the final GOODS LBG sample. 

Postage stamps for each selected galaxy were cut out 
from the GOODS B or V + i images and the SExtrac- 
tor segmentation maps. The pixels assigned to any fore- 
ground objects in the postage stamp were masked out. A 
local sky region was selected to exclude any pixels flagged 
by the SExtractor segmentation map, and the sky value 
within that box was subtracted from the postage stamp. 

3.2. Morphology Measurements 

The elliptical Petrosian radius, Gini coefficient, M20, 
concentration, and average signal-to-noise per pixel 
{{S/N)) were measured from the cleaned, sky-subtracted 
postage stamp of each galaxy. 

A well-defined map of the galaxy's pixels is needed to 
compute G and M20. We choose not to use the isophotal- 
based SExtractor segmentation map. Rather, we recom- 
puted our own map based on the surface brightness of 
the galaxy at the Petrosian radius (rp), where rp is mea- 
sured in elliptical apertures (see LPM04 for a discussion) . 
The elliptical Petrosian radius was calculated assuming 
the SExtractor ELLIPTICITY and THETA values, and 
using the SExtractor computed XCENTER and YCEN- 
TER as a flrst guess. Galaxy pixels with fluxes greater 
than the surface brightness at this radius were assigned 
to the segmentation map of the galaxy. Occasionally this 
new segmentation map is not congruous, and the code 
fails. This generally occurs when either the background 
object masking is not sufflcient or when the object has a 
low surface brightness relative to the sky level (i.e. low 
(S/N)). For the z ~ 1.5 and z ~ 4 samples, any code 
failures were masked by hand and re- analyzed. 



A new estimate of the galaxy's center was then deter- 
mined by minimizing the total second-order moment of 
the pixels within this map: 
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where fi is the flux in each galaxy pixel and Xc, yc is 
the galaxy center. We recalculated r-p using the revised 
center, and a revised segmentation map was used to cal- 
culate the Gini coefflcient and M20. The Gini coefficient 
is defined as the distribution of galaxy pixel fiux values: 
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where n is the total number of pixels in the galaxy seg- 
mentation map. and Xj are the rank-ordered pixel flux 
values (Glasser 1962). The second order moment of the 
brightest 20% of the galaxy pixels is defined as 

M20 = loglO with^ < 0.2 ftot (3) 

Concentration is calculated from the ratio of radii of cir- 
cular apertures containing 80% and 20% of the total fiux, 
where the total flux is the flux within 1.5 circular Pet- 
rosian radii (Conselicc 2003) about the galaxy center de- 
termined by Eqn. 1. Previous definitions of concentra- 
tion have been calculated about the asymmetry center 
(LPM04) or SExtractor center (Conselice 2003). We do 
not calculate asymmetry here due to its instability at 
low (S/N) (LPM04), and instead use the more reliable 
second-order moment center. A bad choice of the center 
can lead to large systematic errors in the concentration 
for very compact objects. 

The average signal-to-noise per pixel {{S/N)) was also 
computed for the pixels in our segmentation map (as 
computed above), using the variance of sky pixel values 
from the previously selected sky region: 
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Because the GOODS ACS images are drizzled and 
remapped to a pixel scale less than the original scale, 
the noise in the sky is correlated. We corrected {S/N) 
for the effects of drizzling using the variance reduction 
factor Fa given in Eqn. A20 of Casertano et al. (2000): 



^correct — 



"sky 

'f7 



(5) 



Fa depends on the ratio s of the output pixel scale to 
the original pixel scale and the drizzled pixel fraction p. 
The GOODS images have s = 0.6 and p = 0.7, giving 
Fa (GOODS) = 0.375. The UDF pixel fraction p is 0, 
therefore Fa (UDF) = 1.0. 

4. SELECTION EFFECTS AND BIASES 

Lyman break galaxies are among the most distant 
galaxies known, and the images we observe (even with 
the HST ACS) are often low signal-to-noise and barely 
resolved. Any analysis of the morphologies of such high- 
redshift objects must also include a careful study of the 
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limits of the morphological measures used. In LPM04, 
wc tested the stability of G, M20, and C with decreasing 
(S/N) and spatial resolution for local galaxies observed 
in rest-frame R. First, the images for set of 7 galaxies 
of different morphological types were held a fixed spatial 
resolution (better than 120 pc per resolution element) 
but added to increasingly higher Poison noise images to 
measure the effect of noise on the morphological quan- 
tities. Then the images were smoothed to decreasing 
spatial resolutions but scaled such that (S/N) was con- 
stant to measure the effect of low spatial resolution on 
the morphological quantities. These tests showed that for 
relatively nearby galaxies observed in rest-frame optical 
wavelengths, G and M20 can be recovered for observa- 
tions with (S/N) > 2 — 3 and spatial resolutions better 
than 1 kpc. 

However, other effects may come into play when mea- 
suring the morphologies of distant galaxies. The local 
sample of LPM04 may not adequately reflect the range 
of morphologies or physical scales of objects at 2 > 1. 
High redshift galaxies may be appear more irregular, 
due to both evolution and morphological K-corrections 
(Abraham et al. 1996). Also, high redshift galaxies are 
physically smaller than local galaxies (e.g. Ferguson et 
al 2004). Finally, the morphologies of galaxies smaller 
than a few times the image's PSF cannot be measured, 
because most of the object's spatial information is lost. 

In this work, we use three approaches to to determine 
the morphological selection effects and biases of our dis- 
tant galaxy samples: 1) a comparison of the measured 
morphologies for large sample of galaxies observed by 
both the GOODS survey and the deeper UDF survey; 

2) a series of simulations of face-on smooth exponential 
disks and bulges with varying sizes and magnitude; and 

3) simulations in which we artificially redshift our 2; ~ 1.5 
sample to 2; = 4. The first test will help us understand 
at what (S/N) level the measured morphologies of real 
(and morphologically complex) distant galaxies are re- 
producible. The second test will help determine any off- 
sets from the "true" morphologies as a function of size, 
surface brightness, and light profile. The last test will 
allow us to directly compare the observed morphological 
distributions at z ^ 1.5 and z ^ 4. 

4.1. GOODS V. UDF morphologies 

One of the advantages of the GOODS data-set is that 
a sub-sample of the GOODS CDFS field has been re- 
observed with HST ACS in the same bandpasses to a 
substantially greater depth for the Ultra Deep Field Sur- 
vey. By comparing the size, magnitudes, and morpholo- 
gies calculated from the extremely high signal-to-noise 
images from the Ultra Deep Field survey to those calcu- 
lated for the exact same objects in the shallower GOODS 
images, we can learn a great deal about the reliability of 
the morphologies measured for the entire GOODS sam- 
ple. 

We have selected ^ 750 objects from the UDF catalog 
with z magnitudes brighter than z{FSbQW) = 28 mag 
as our control sample. No color or redshift cut was ap- 
plied to this sample, and it includes galaxies as bright 
as z{FSbQW) = 20 mag. Both the UDF and overlap- 
ping GOODS summed V + i images were pre-processed 
as described above without modification to the SExtrac- 
tor segmentation maps, and the morphologies and {S/N) 



for all 750 objects were measured in both the GOODS 
and UDF images. From this sample, we selected galax- 
ies with UDF (S/N) > 3.0 and UDF elliptical Petrosian 
radii Vp > 0.3 ". We also removed objects for which the 
morphology code failed, i.e. objects for the Gini segmen- 
tation map was not congruous. The morphology code 
failed for less than 5% of galaxies brighter than V + i = 
26.5 in the GOODS-S images. 

The differences between the GOODS and UDF mor- 
phology and size measurements for the resulting 571 ob- 
jects are plotted as a function of the GOODS (S/N) 
in Figure 1. The median offset (GOODS-UDF) and 
standard deviations for the 67 objects with 2.5 < 
(S/AT) (GOODS) < 3.5 is 

Arp = -0.057" ± 0.097" 

AG = -0.018 ± 0.033 , . 

AM2o= +0.062 ±0.100 

AC = -0.219 ± 0.448 

We repeated this analysis for the B{F435W) GOODS 
and UDF images and found very similar offsets and stan- 
dard deviations at 2.5 < (S'/A') (GOODS) < 3.5. At 
(S'/A^) (GOODS) < 2.5, the differences between the mea- 
sured UDF and GOODS morphologies show large scat- 
ter. Therefore, we will restrict our analysis of galaxy 
morphologies to objects with {S/N) > 2.5 

At low (S/N), the Petrosian radii measured in the 
GOODS images arc; slightly smaller than those measured 
in the UDF imagers. This is likely to be the result of light 
in the outer low surface brightness regions being lost in 
the sky noise. This bias toward smaller sizes is probably 
the reason for a small offset in G at low {S/N). The 
lowest surface brightness pixels are not included in the 
segmentation map, and the Gini coefficient is artificially 
lowered. The same bias may also be the cause for the 
mean offset to higher M20 values at low {S/N). M20 
is the ratio of the second order moment of light for the 
bright regions of the galaxy to the total second order mo- 
ment for the entire galaxy. This total moment depends 
on the Petrosian radius; if the Petrosian radius is artifi- 
cially lower, then the total moment decreases while the 
moment of the brightest region remains unaffected. The 
concentrations measured in the GOODS images appear 
to be systematically lower than those measured in the 
UDF images, even at relatively high {S/N). This offset 
(—0.22) is small relative to the dispersion (0.45) and the 
range of possible concentrations. This is also the result of 
missed light in the outer regions of the galaxy - we mea- 
sure systematically smaller apertures containing 80% of 
the total light in the GOODS images, even for galaxies 
with {S/N) ~ 4. 

4.2. Bulge and Disk Simulations 

To determine the robustness of the measured 
morphologies at smaller sizes and fainter fiuxes, 
we have simulated a series of face-on exponential 
disks and r^/'* bulges of varying sizes and magni- 
tudes. The simulated galaxies are constructed using 
IRAF.ARTDATA.MKOBJECTS task and are added to 
a blank region of the V + i or B images. We detect and 
measure the simulated objects in the same way as for 
our real galaxy samples. A noise-free image of de Vau- 
couleurs bulge with a half-light radius of 600 resolution 
elements has G = 0.600, M20 = -2.47, and G = 4.34, 
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while a noise-free image of an exponential disk with the 
same half-hght radius has G = 0.473, M20 = —1.80, and 
C = 2.71 (horizontal dotted lines in Figure 2). 

We use these simulations to estimate any 
size/luminosity dependent biases in the morphol- 
ogy measurements. At z > 3, the ACS images have 
physical scales ~ 600 pc per PSF FWHM (0.12 "). 
Therefore, we expect that intrinsically small and/or 
highly concentrated galaxies for which the central 
regions are not resolved will have systematic biases in 
their measured morphologies. Wc find that C, M20, and 
to a lesser extent G are affected. In Figure 2, we plot 
the measured morphologies and sizes for both bulges 
and disks detected with {S/N) > 2.5, with different 
grey-scales indicating different total input magnitudes. 
Our simulations show that disk-dominated galaxies with 
rp < 0.3" arc not sufficiently resolved to give accurate 
G. At small angular sizes, the measured G approaches 
that expected for a point source (G ~ 0.7). We find 
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that M20 is biased low for the bulges with sizes less 
than 1 "; local well-resolved spheroids have M20 ~ — 2.5 
(LPM04), while our simulations produce M20 ~ —2.0. 
Nevertheless, bulges and disks still have distinguishable 
M2a values. At large sizes and faint magnitudes, we 
find that the measured radii are smaller than their true 
radii as the flux from the outer low surface brightness 
regions is lost. This is a problem particularly for bulges. 
Because the lowest surface brightness pixels are lost, 
this effect can also artificially lower the measured G by 
0.05-0.06. 

Concentration is unreliable at small sizes, because the 
inner aperture containing 20% of the total light not re- 
solved. The observed step-like behavior of concentration 
for the simulated bulges at Vp < 0.7" is a result of inter- 
polating the inner 20% light radius (^20) at less than 1 or 
2 pixels. At radii < 1 pixel, the interpolated r2o is likely 
to be an overestimate because of the steep inner surface 
brightness profile of r^/* bulges, and results in an under- 
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Fig. 2. — The measured morphologies as a function of measured size (rp) and magnitude for simulated face-on exponential disks (circles) 
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horizontal dotted lines show the "true" morphological values. The dashed line show the rp = 0.3" cutoff. 



estimate of C. The estimated r2o approaches the true 
value at r2o = 1 pixel, and then is again overestimated 
at r2o > 1 pixel. Once the measured r2o is greater than 
2 pixels, this effect is not as significant. 

In addition to selecting galaxies with {S/N) > 2.5, 
we add the additional selection criterion that the galaxy 
must have measured rp > 0.3 ", because smaller objects 
are likely to have artificially high G values. This will also 
effectively eliminate any possible stellar contaminants in 
our LBG sample. The net result is that our sample will 
span a limited range in sizes, as larger objects are likely 
to have too low {S/N) and smaller objects are too com- 
pact to be measured. In Figure 3, we show the measured 
rp and isophotal magnitudes for the GOODS and UDF 
B-dropouts and the GOODS z ^ 1.5 sample as a func- 
tion of (S/N). All panels in Figure 3 are plotted on 
the same physical scale. The high {S/N) z ~ 1.5 galax- 



ies are intrinsically larger and fainter than the GOODS 
B-dropouts, but overlap in absolute magnitude and size 
with the UDF B-dropouts. The redshift distribution of 
the selected z 1.5 subsamples is shown in Fig. 4. 

4.3. Artificial redshift simulations 

In Figure 3, we demonstrated that the UDF z ^ 4: 
and GOODS 2 ~ 1.5 samples span a similar range of 
intrinsic sizes and rest-frame FUV luminosities. There- 
fore, any differences between the morphological distri- 
butions of these two samples is likely to be the re- 
sults of intrinsic differences in the samples as opposed 
to selection effects. However, the selection functions are 
not identical. Therefore we have selected those objects 
from the z ~ 1.5 sample with measurable morphologies 
{rp > 0.3, {S/N) > 2.5) and artificially redshifted these 
to z; ~ 4. The redshifted images were given the same ex- 
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isophotal magnitude for GOODS and UDF B-dropouts (upper panels) z ~ 1.5 sample and artificially redshifted 



z ~ 1.5 — » 4 sample (lower panels). All four panels are set to the same physical scale in absolute magnitude and intrinsic size (kpc). 



posure times and sky background as the UDF V + i im- 
ages. Wc find that z 1.5 objects with isophotal surface 
brightnesses < 23.6 mag per sq arcsec in the GOODS 
B images drop below the UDF V + i limiting surface 
brightness of ~ 26.1 mag per sq arcsec when artificially 
redshifted because of the (1 + z)^ surface brightness dim- 
ming. Therefore only 23 of the artificially redshifted 
objects have {S/N) > 2.5. For these objects, the off- 
sets in morphology between the original and redshifted 
images are similar to those found in Equation 6. The 
measured Vp and isophotal magnitudes for the artificially 
redshifted sample is shown in Figure 3. The UDF z ^ A 
sample probes a volume roughly four times that of the 
un-redshifted GOODS z ~ 1.5 sample, assuming a limit- 
ing magnitude of ^ + i ~ 27 and B ^ 25 for each sample 
respectively. 

5. GALAXY MORPHOLOGIES AT Z ~ 1.5 AND Z ~ 4 



In Figures 5-7, we show the G, M20, and C values for all 
galaxies with (S/N) > 2.5 and rp > 0.3"(small crosses), 
as well as the artificially redshifted sample. We find that 
the 2 ~ 4 LBGs have median G ^ 0.58, M20 ^ -1-6 
and C ~ 3.8. We note that Ferguson et al. (2004) 
found similar concentration values for the ACS i;-band 
observations of GOOD LBGs, but that Conselice et al. 
2004 and LPM04 reported significantly lower concentra- 
tion values for the observed near-infrared morphologies 
of 2 < z < 3.5 galaxies in the Hubble Deep Field North. 
We suspect that the larger PSF of the NICMOS images 
has biased the measured concentrations to lower values. 
The mean G and M20 values are very similar to those 
found by LPM04 for the rest-frame u and B morpholo- 
gies of z ~ 2 and 2: ~ 3 galaxies in the Hubble Deep Field 
North. 

A series of two-dimensional Kolmogorov-Smirnov (K- 
S) tests (Fasano & Franceschini 1987) give a less than 
37% probability that bright GOODS LBGs and the 
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Fig. 1. Spectroscopic rcdshift histogram for GOODS 1.2 < 
z < 1.8 sample (shaded histogram; 55 galaxies with {S/N) > 2.5, 
rp > 0.3" selected for the morphological analysis) . 



fainter UDF LBGs arc drawn from the different distri- 
butions in G, M20 and C. The z ~ 1.5 sample's median 
morphological values are slightly different, with median 
G ~ 0.55 (with a tail to lower G), M20 ^ —1.5, and 
C ~ 3.3. These median values do not change when the 
z ^ 1.5 sample is artificially rcdshiftcd to z = 4. The 
K-S statistical probability that the morphological distri- 
butions of the artificially-redshifted ^ ~ 1.5 ^ 4 sample 
is different from the UDF .2 ~ 4 LBG sample is greater 
than 91.5%. 

The postage stamps for a subset of galaxies in both rcd- 
shift ranges are shown in Figures 8-10. In all samples, 
galaxies with M20 > — 1.1 have well-separated double or 
multiple bright nuclei. Galaxies with M20 < —1.6 are 
relatively smooth with a single nucleus. Objects with 
intermediate M20 values can be irregular in appearance, 
often with a bright nucleus with tidal tail or fainter knots. 
Galaxies with G > 0.6 have very bright nuclei, and in- 
clude both smooth "bulgc"-like objects and close pairs. 
In the lower-redshift sample, we observe a number of 
galaxies with G < 0.55, M20 > -1-5, C < 3 . These tend 
to have extended, lower-surface brightness star-forming 
regions, reminiscent of local disk galaxies (Figure 10). 

5.1. FUV Morphological Classification 

Although we are able to robustly measure the rest- 
frame FUV morphologies of our z ^ 1.5 and ^ 4 sam- 
ples, we cannot simply apply the quantitative classi- 
fication criteria that was originally developed for the 
optical morphologies of local galaxies. Young star- 
forming populations completely dominate the rest-frame 
FUV morphologies of local galaxies, while old or dust- 
enshrouded stellar populations all but disappear. Thus 
the FUV morphologies of normal galaxies appear sig- 
nificantly more disturbed than their optical morpholo- 
gies (Giavalisco et al. 1996, Hibbard and Vacca 1997; 
Goldader et al. 2002) . In order to re-calibrate our merger 
classification criteria, we have examined the FUV mor- 
phologies of four bright well-resolved galaxies observed 



by the GALEX Nearby Galaxy Survey. While a detailed 
study of a much larger sample of galaxies is needed, our 
preliminary work here in combination with our bulge and 
disk simulations provides a check on our revised FUV 
morphology classifications. 

The GALEX FUV images of NGC1399 (El), NGC1068 
(Sb), NGC2403 (Scd), and NGC520 (Arp 157, merger) 
have been artificially redshifted to z = 1.5 and z = 4, 
convolved with the ACS PSF simulated using Tiny Tim, 
and added to the GOODS B, V + i, and UDF V + i 
sky background (Figure 8). Because the absolute Mpuv 
of the local galaxies are typically > —16 and would be 
undetectable at z > 1.5, the artificially redshifted galax- 
ies have been brightened to Mpuv = —21. We have 
not applied any K-corrections to the images. We have 
measured the morphologies of the artificially redshifted 
galaxies as they appear at z = 1.5 and z = 4 and plotted 
their G, M20, and C values in Figure 5-7. The elliptical 
NGC1399 (circle) has G > 0.6, M20 < -2.0 and C > 4.3, 
in rough agreement with our r^^^ bulge simulations (up- 
per right-hand box). The Sb NGC1068 (triangle) has 
G ^ 0.58, M2Q ^ -2.0, and C > 4.3, which also agrees 
with our bulge simulations. The Scd NGC2403 (square) 

has G ~ 0.55, M20 1.4, and C ~ 2.9. This gives 

a slightly higher G and M20 values than our smooth 
exponential disk simulations (lower left-hand box), as 
might be expected for a real spiral with bright star- 
forming knots. The merger NGC520 (star) has G ~ 0.59, 

M20 1.0, and C ~ 4.3. NGC520 has been classified 

as an intermediate stage merger of two gas-rich spirals 
(Hibbard & van Gorkum 1996), and only one of its nu- 
clei is visible in the FUV. 

Based on the redshifted GALEX images, our bulge and 
disk simulations, and visual inspection, we have clas- 
sified the z ~ 1.5 and z ~ 4 FUV morphologies into 
three types. Major-merger candidates such as close pairs 
and/or objects with multiple bright nuclei or bright tidal 
tails have M20 > — 1.1 (left shaded region of Figure 5). 
Bulge-dominated systems (E-Sb) are those with G > 0.57 
and M20 < — 1.7 (upper right shaded region of Figure 5). 
Galaxies which do not match any of these criteria are 
called "transition" objects. We stress that this initial 
classification scheme is more conservative than that of 
LPM04 because of the lack of a large calibration sam- 
ple of FUV galaxy images and will only identify major 
mergers in their most morphologically-disturbed stage. 
We expect that some of the transition objects will be 
minor-mergers and post mergers, as well as late-type 
star-forming disks which appear knotty in the rest-frame 
FUV. We find that C-based classifications often disagree 
with those based on G — M20, as C is less robust at small 
sizes and is less sensitive to substructure (Figures 6 and 
7). 

5.2. Merger and Spheroid Fractions at z ^ 1.5 and 

In Table 1 , we give the observed merger and bulge frac- 
tions for our UDF z - 4 and GOODS z - 4 LBGs, the 
GOODS z ^ 1.5 star-forming galaxies, and the artifi- 
cially redshifted sample. At Mpuv < —20, the observed 
fraction of merger candidates is similar for all samples (~ 
10%: - 20%). In Figure 3, we showed that the GOODS 
and the UDF LBG samples probe different parts of the 
luminosity function, while the UDF and z ~ 1.5 samples 
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are better matched. In this section, we wih use the com- 
bined GOODS and UDF LBG samples to constrain the 
morphological distribution of 2: 4 star-forming galax- 
ies, but use only the UDF z ^ 4 and z ~ 1.5 samples to 
constrain any evolution with redshift. 

If we classify all objects with M20 < —1.1 as on-going 
major merger candidates, the observed major merger 
fraction at z ~ 4 is ~ 14% for GOODS Bdrops (5/36) 
and ~ 10% for UDF Bdrops (5/46). However, we have 
not corrected this fraction for incompleteness. In order 
to do a rough incompleteness correction, we have visually 
inspected the LBGs with (S/N) < 2.5 and V + i magni- 
tudes brighter than 25. These objects have rp > 0.5"and 
were excluded from our original sample because of their 
low surface brightness and {S/N). Fifteen out of the 
25 bright excluded objects were visually classified as dis- 



turbed, therefore it is possible that we are missing a sig- 
nificant fraction of large, bright but lower-surface bright- 
ness merger candidates. Four of these objects were visu- 
ally classified as bulge- dominated. We have also excluded 
a large number of low-surface brightness GOODS LBGs 
at 25 < V + j < 26. Examination of the UDF LBGs 
at the same magnitude and sizes indicates that ~ 22% 
of the excluded GOODS LBGs with 25 < F -|- z < 26 
are likely to be major merger candidates and ^ 22% are 
likely to be btdge-dominated. If we correct the GOODS 
LBG merger fraction for these rough incompleteness esti- 
mates, we find that average fraction of major- merger can- 
didates is ~ 25% for all LBGs with V + i < 26, rp> 0.3". 
If we consider only the brightest LBGs with V + i < 25, 
rp > 0.3", the major-merger fraction could be as high as 
~ 50% once incompleteness effects are included. How- 
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ever, our corrections are highly uncertain, we do not ac- 
count for LBGs which are undetected, and the number 
of bright LBG is small (16). Nevertheless, these may be 
considered rough upper limits to the morphologically- 
determined merger fraction for star-forming galaxies. 

It is unclear whether our major-merger fraction for 
2; ^ 4 LBGs is consistent with previous high-redshift 
quantitative morphology studies based on the rotational 
asymmetry parameter. Consclicc ct al. (2003) found a 
50% merger fraction for Mb < —21 galaxies at 2: ~ 2 — 3 
based on a sample of 10 objects in the HDFN. If we 
assume {FUV — B) ^ 1, then our merger fraction es- 
timate of 25% is about a factor of two lower for the 
same magnitude range [Mpuv < —20) based on a 
sample of 60 LBGs in GOODS and UDF. However, at 
Mpuv < —21 {Mb < —22) our visual classification sug- 
gests that ~ 50% are disturbed objects, so there may 
be a trend towards higher merger fractions at brighter 
luminosities. 



We find that a significant fraction of 2; ^ 4 objects 
have morphologies consistent with spheroidal morpholo- 
gies: 11/36 of the GOODS LBGs and 19/46 of the UDF 
LBGs lie within the spheroid criteria in G — M20 space. 
Brighter than Mpuv ~ -20, roughly 30% of the GOODS 
and UDF galaxies appear to be bulges. Many of the UDF 
spheroids are fainter than Mpuv = —20 (11/22). This 
is in contrast to the Mpuv > —20 z ~ 1.5 galaxies, of 
which only 2/15 of the artificially redshifted sample are 
bulge-Hke. 

The remaining ~ 50% of the galaxies in the high- 
redshift samples have higher M20 and G than expected 
for smooth disks or bulges but lower M20 values than 
the major-merger candidates. Some of these objects lie 
close in G — M20 to the artificially-rcdshifted FUV im- 
age of late-type spiral NGC2403, therefore many may 
be star- forming disks. In the local universe at optical 
wavelengths, higher M20 and G values are indicative of 
merger signatures such as very bright nuclei and tidal 
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tails, and it is also possible that a fraction of the transi- 
tional objects are minor-mergers or post-merger systems 
which are disturbed but do not possess bright multiple 
nuclei and high M20 values. 

The z ^ 1.5 objects are more likely to fall into this 
transitional class than the UDF LBGs (60% of the z ^ 
1.5 sample and 70% of the artificially redshift sample). 
Also, the G distributions of the transitional objects are 
slightly different at the two epochs: less than 5% (1/46) 
of the observed UDF z ~ 4 sample has G < 0.5, as 
opposed ~ 20% at lower redshift ( 10/55 of the ^; ~ 1.5 
sample and 5/23 of the artificially redshifted 2: ^ 1.5 ^ 4 
sample). Low G values usually indicate disk-like galax- 
ies with more uniform surface brightness profiles. Is 
the lack of low G galaxies at high redshift the result of 
an observational bias? We clearly select against large 
disks because of their low-surface brightnesses. How- 
ever the lower-redshift sample has low G objects with 
sizes, total magnitudes, and effective surface brightnesses 



which should have been detectable in the UDF (Fig 3). 

Our artificial redshift simulations suggest that the de- 
tected fraction of low G objects should not change due 
to redshift-dcpendcnt selection effects. 

The slightly lower Gini coefficients, concentrations, 
and second-order moments at 2: ~ 1.5 may be partially 
a result of larger dust obscuration. Dust may have sig- 
nificant effect on the observed morphologies in the FUV 
than at longer wavelengths. In a study of local ULIRGs 
with HST STIS images in the FUV/NUV and WFPC2 
and NICMOS images at optical and near-infrared wave- 
length, Goldader et al. (2002) found than many star- 
forming knots were obscured in the FUV/NUV, result- 
ing in substantial changes in the observed morphologies 
and light distributions. The z ~ 4 Lyman break galaxies 
may not suffer from large dust extinction because Ly- 
man break technique automatically selects blue objects 
with low to moderate extinctions; however, the 2; ~ 1.5 
sample is likely to possess dust and have higher internal 
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color dispersions. In a study of the internal FUV-optical 
color dispersions of Hubble Deep Field North 2 ~ 2.3 
and 2 ~ 1 samples, Papovich et al. (2005) found that 
the lower rcdshift sample had significantly more internal 
color dispersion than the high redshift sample, a result 
of either patchier dust distributions or larger spatial dis- 
persions in the age of the stellar populations. 

5.3. FUV - B Colors and UV Slope 

The southern GOODS field (CDFS) has deep ground- 
based near-infrared imaging in J, H, Kg with the ISAAC 

camera at the VLT (Vandame et al 2004). The GOODS 
ACS images were degraded to the ISAAC image reso- 
lution, and matched aperture photometry was obtained 
for the CDFS/UDF galaxies (Dahlen et al. 2005). Eight 
GOODS-S LBGs and twelve UDF LBGs (one of which is 
also in the ISAAC-detectcd GOODS LBG sample) from 
our (S/N) > 2.5, Vp > 0.3 " were detected in the ISAAC 
Ks image and 16 GOODS-S 2; ~ 1.5 galaxies were de- 



tected in the ISAAC J image. We have examined the 
rest frame B isophotal magnitudes {K and J for .2 ~ 4 
and 2; ~ 1.5 samples), rest-frame FUV — B isophotal 
colors, and rest-frame FUV — NUV isophotal colors as 
a function of galaxy morphology in Fig. 12 and 13. At 
z ^ 4, the z - K color corresponds to ~ 1700A - 4400A, 
while at 2 ~ 1.5, the B — J color corresponds to ~ 1740A 
- 44OOA. We use the rest-frame FUV — NUV colors as a 
proxy for the UV spectral slope /3 (Fx oc A'', Calzetti et 
al. 1994), using the ISAAC-degraded z — J isophotal col- 
ors (~ rest-frame 1700 A - 2200A) for the 2^4 galaxies 
and using the ACS B — V isophotal colors (~ rest-frame 
1740A - 2500 A) for the z ~ 1.5 galaxies. 

We find no strong correlation of rest-frame 4400A lu- 
minosity with FUV galaxy morphology for either sample 
(upper panels, Fig. 12 and 13). We find that majority of 
z ^ 4 LBGs have z — K < 1.5. There are two red LBGs 
with i — K > 2.5; such red rest-frame FUV — B colors 
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Fig. 9. — V + i 1.5" x 1.5" images for a subset of GOODS 2: ~ 4 LBGs. Classifications are based on Figure 5. 
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Fig. 10. — V + i 1.5" X 1.5" images for a subset of UDF 2; ~ 4 LBGs. Classifications are based on Figure 5. 



indicate strong 4000A breaks and relatively old stellar 
populations (> 0.5 — 1 Gyr), while the bluer colors of 
the majority of the z ^ 4 LBGs indicate younger stellar 
populations. The blue LBGs show no obvious correla- 
tion of rest-frame FUV — B color or UV slope (z — J or 
rest-frame FUV - NUV) with G, M20 or C. The two 
red LBGs with i — k > 2.0 have morphologies consistent 
with spheroids (G > 0.5, M20 < —1.5) and arc also the 
brightest objects in the iiT-band sample. The z ~ 1.5 
sample has similar FUV — B colors but show steeper 
UV slopes {B — Vov rest-frame FUV — NUV), implying 
either larger dust extinctions or older ages. There are no 
obvious color-morphology trends. The reddest object in 



the 2; ~ 1.5 sample is also consistent with a bulge-like 
morphology (G = 0.60, M20 = -1.7, G = 5.1). 

6. DISCUSSION 

We find no strong evolution in the rest-frame FUV 

morphologies of star-forming galaxies between z ^ 4 and 
z 1.5, with a difference in look-back times of roughly 
3 billion years. Wc see subtle changes in the fraction of 
spheroids and low G disk-like objects, but the fraction of 
objects with bright multiple nuclei is indistinguishable 
for the two samples. At first glance, this result is some- 
what surprising given the large time difference between 
the two observed epochs and the apparent increase in 
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Fig. 11. — B 1.5" X 1.5" images for a subset of GOODS z ~ 1.5 galaxies. Classifications are based on Figure 5. 



TABLE 1 



Morphology Distributions 






Mpuv 


< -21 


-21 < Mpuv < -20 


-20 < Mpuv < -19 


Sample 


Type 


Nobs 


Fraction 


Nobs 


Fraction 


Nobs 


Fraction 


GOODS z ~ 4 


Total 


11 




25 










Mergers 


2 


0.47 (0.18) 


3 


0.21 (0.12) 








Bulges 


3 


0.19 (0.27) 


8 


0.23 (0.32) 








Other 


6 


0.34 (0.55) 


14 


0.56 (0.56) 








UDF 2; ~ 4 


Total 


5 




19 




22 






Mergers 


3 


(0.6) 


1 


(0.05) 


1 


(0.05) 




Bulges 


1 


(0.2) 


7 


(0.37) 


11 


(0.50) 




Other 


1 


(0.2) 


11 


(0.58) 


10 


(0.45) 



GOODS z ~ 1.5 Total 3 • • • 23 ■ • • 28 

Mergers 1 (0.33) 4 (0.17) 5 (0.18) 

Bulges 2 (0.66) 4 (0.17) 8 (0.30) 

Other (0.00) 15 (0.65) 16 (0.57) 





GOODS 2 


~ 1.5 Total 


2 




8 


13 




Mergers 


1 


(0.5) 


1 (0.125) 


(0.0) 




Bulges 





(0.0) 


1 (0.125) 


2 (0.15) 




Other 


1 


(0.5) 


6 (0.75) 


11 (0.85) 


Note. — 


The fractions in parentheses 


are 


the number of g 


alaxies in the each morphological bin 


divided by the total measured objects in each 



sample and have not been corrected for incompleteness. 



dust extinction at later times. 

However, the dominant mode of star formation at both 
epochs is expected to be cohisional starbursts triggered 
minor-mergers - 45% of p* at 2: ~ 1.5 and 70% of at 
2 ~ 4 (SPFOl). Therefore many FUV-bright galaxies at 
both epochs arc likely to be minor mergers which exhibit 
a dominant bright nucleus with a disturbed morphology. 
In fact, this is what we find for both our z ~ 1.5 and 
z ^ 4 samples. Over half of the galaxies have been clas- 



sified as "transition" objects between the spheroidal mor- 
phologies and the double-nucleated major-merger candi- 
dates. While it is likely that some of these transition 
galaxies are post-major mergers or knotty star-forming 
disks, a significant fraction are likely to be minor merg- 
ers. 

We also find very similar fractions of major merger can- 
didates at both epochs. Our UDF z ~ 4 major-major 
candidate fraction of observed star-forming galaxies is 
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Fic;. 12. Isophotal K AB magnitudes, z — K colors (~ rest-frame 1700A - 4400A), and z — J colors (~ rest-frame 1700A - 2200A) v. 
morphology for GOODS-South (crosses) and UDF (circles) ^ ~ 4 samples. 



~ 10% for Mpuv < —19, while the 2; ~ 1.5 major-merger 
candidate fraction is ~ 10 — 20% over the same luminos- 
ity range. Wc caution that these observed fractions are 
not "merger rates" ; in order to calculate a true merger 
rate, we must correct for the selection biases in our sam- 
ple against low-surface brightness or low star-formation 
objects as well as the probability and timescales for ob- 
serving mergers with M20 < — 1.1. As we discussed in 
§2, identifying non-star-forming galaxies at these epoch 
is challenging, therefore we limit our analysis to star- 
forming galaxies only. The probability of observing a 
merging system during a period of disturbed morphol- 
ogy will depend on a number of factors including the 



initial conditions of the merger system (mass ratio, gas 
fraction, orbital parameters), the viewing angle, and the 
timescale of the merger. A large suite of hydro-dynamical 
merger simulations required to calculate this probability 
are currently being analyzed (Cox et al. 2005, Jonsson ct 
al. 2005, Lotz et al. in prep). From our initial analysis of 
gas-rich major merger simulations, we expect equal mass 
merging spirals on a prograde orbit with a gas-fraction 
equal to 0.52 to have M20 > -1-1 for - 500 Myr. 

The number of major mergers per halo per Gyr is pre- 
dicted to be ~ 15% at 2 ~ 1.5 and ~ 25% at 2 ~ 4, 
where the minimum halo circular velocities are > 50 km 
s~^ (Gottlober, Klypin, & Kravtsov 2001). If number 
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of detected star-forming galaxies per typical halo is 2 
and does not change strongly from 2~4to2;~1.5, then 
our results are roughly consistent with the predicted ma- 
jor merger rate for halos. Given that our samples are 
fairly small, and that we have not adjusted the observed 
major merger fractions for undetected non-star- forming 
galaxies, it is surprising that our numbers are so con- 
sistent with theory. Clearly a more detailed comparison 
between the merger rates and morphology distribution 
between hierarchical models and the data is needed that 
takes into account the many selection effects and biases. 

At later times, quiescent star-formation should become 
increasingly important. SPFOl predict that 40% of 



star-formation density occurs in a quiescent (disk-star 
formation) mode at 2 ~ 1.5 as compared to ~ 20% at 
z ^ 4. This suggests that the fraction of FUV-bright 
disk galaxies will increase at later epochs. While we find 
relatively few "smooth" disks at cither epoch, there are 
more objects with G < 0.5 at 2; ~ 1.5 which are con- 
sistent with low-surface brightness disks (see Figure 11) 
than at z 4. These disk candidates make up ^ 20% of 
the entire z ~ 1.5 sample and a third of the transition 
objects at that epoch. 

Finally, the observed fraction of star-forming galaxies 
with spheroidal morphologies is similar at both epochs. 
However, when the z ~ 1.5 sample is artificially red- 
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shifted to z ^ 4 the fraction of spheroids in the faintest 
bin is significantly lower than that observed in the UDF 
2 ~ 4 sample. If these bulges are only formed as the end- 
stage of a major merger, then similar fractions would be 
expected. However, if minor-mergers and interactions 
also produce blue bulges (Kannappan, Jansen, & Barton 
2004), then one might expect significantly higher rates 
of star-forming spheroids at higher redshift given the in- 
creased importance of minor- mergers at those times. On 
the other hand, there may not be enough time before 
2 ~ 4 for significant numbers of proto-bulges to form via 
minor-mergers, as the timescale for a minor merger is > 
2 Gyr (Cox 2004). 

Thus the observed distributions of star-forming galaxy 
FUV morphologies and the weak evolution between z ^ 
1.5 and 2; ~ 4 appears to be generally consistent with 
rates of major and minor mergers and the importance 
of quiescent star-formation predicted at those epochs by 
hierarchical galaxy evolution models. However, we cau- 
tion that we have largely ignored the morphological ef- 
fects of dust and the evolution of typical dust extinc- 
tion with redshift. The FUV — NUV colors suggest 
that the z ~ 1.5 sample is more extincted than the the 
LBGs, as might be expected if the galaxies are chem- 
ically evolving with time. Strong dust extinctions are 
likely to hide regions of intense star formation, and shift 
observed morphologies to lower G and possibly lower C 
and M20 values, as well as remove some heavily-extincted 
objects from our samples altogether. Nevertheless, since 
the observed morphologies does change significantly from 
2 ~ 4 and z ~ 1.5, increased dust does not appear 
to dominate the observed FUV morphologies for mod- 
erate extinctions. We also caution that our conclusions 
here only apply to galaxies which have sufficient star- 
formation to meet the LEG color-criteria at 2; ~ 4 and 
to have emission-line determined redshifts at z 1.5. 
We are most certainly missing galaxies at both epochs 
with no or strongly enshrouded star-formation. 

7. SUMMARY 

1) We find no significant difference between the rest- 
frame FUV morphological distributions of GOODS and 
UDF 2; ~ 4 LBGs samples. Lyman break galaxies with 
morphologies consistent with on-going major mergers 
make up ~ 10 - 25% of the Mpuv < -20 GOODS 
and UDF selected samples. A significant fraction of 
the LBGs appear to have relatively undisturbed bulge- 
like FUV morphologies (~ 30%). The remaining ^ 
50% have G and M20 values higher than expected for 
smooth exponential disks or bulges but do not possess 



double nuclei, and could be a mixture of star-forming 
disk, minor-mergers, and post-merger systems. The ob- 
served major-merger fraction and large number of minor- 
merger candidates are consistent with current numerical 
(Gottlober et al. 2001) and semi-analytic model predic- 
tions (Somerville et al. 2001) for high redshift galaxies. 

2) The star-forming z ^ 1.5 sample has rest-frame 
FUV morphologies which are not significantly different 
from the UDF 2: ~ 4 galaxies, with similar fractions 
of major-merger candidates. The lower redshift sam- 
ple does possess slightly more objects with lower G and 
slightly lower C and higher M20 values, implying more 
extended, lower surface brightness disk-like star-forming 
regions, as well as a smaller fraction of faint spheroids. 
This is also roughly consistent with the major-merger 
rates and the increased importance of quiescent (disk) 
star-formation at 2; ~ 1.5 predicted by current hierarchi- 
cal models. 

3) The two LBGs with rest-frame FUV - B colors 

consistent with strong 4000A breaks and > 500 Myr stel- 
lar populations have bulge-like morphologies, as does the 
reddest 2^1.5 galaxy. The UV slopes of the lower red- 
shift sample are significantly steeper than those of the 
2 ~ 4 LBGs. Although dust extinction may be more 
important in the z ~ 1.5 samples, the apparent lack of 
evolution in FUV morphology since 2 ~ 4 suggests that 
moderate extinctions do not strongly bias the observed 
FUV G and M20 values. High-resolution rest-frame op- 
tical morphologies for large samples of high redshift ob- 
jects are needed to constrain the morphologies of galaxies 
which are heavily dust-extincted or not forming stars. 
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